Photodissociation and radiative association of HeH + in the metastable triplet state 
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We investigate the photodissociation of the metastable triplet state of HeH + as well as its forma- 
tion through the inverse process, radiative association. In models of astrophysical plasmas, HeH + 
is assumed to be present only in the ground state, and the influence of the triplet state has not 
been explored. It may be formed by radiative association during collisions between a proton and 
metastable helium, which are present in significant concentrations in nebulae. The triplet state can 
also be formed by association of He + and H, although this process is less likely to occur. We compute 
the cross sections and rate coefficients corresponding to the photodissociation of the triplet state by 
UV photons from a central star using a wave packet method. We show that the photodissociation 
cross sections depend strongly on the initial vibrational state and that the effects of excited elec- 
tronic states and non-adiabatic couplings cannot be neglected. We then calculate the cross section 
and rate coefficient for the radiative association of HeH + in the metastable triplet state. 



I. INTRODUCTION 

The helium hydride ion HeH + is one of the most ele- 
mentary molecular ions and the first to form in the early 
universe [l[ , through the direct radiative association pro- 
cess 0] 

He(ls 2 )+H+ — > HcH+(X X E + ) + hv , (1) 

which occurs in the ground state of the molecular ion. 
In addition, HeH + can also be formed by a spontaneous 
radiative transition from the vibrational continuum of 
the excited A : S + state to a discrete vibrational level of 
the ground state [1, 0] : 

He+(ls) + H(ls) — ► HcH+pf + hv (2) 

These mechanisms of production of HcH + were inves- 
tigated along with destruction processes in order to esti- 
mate the abundance of the molecular ion in the early uni- 
verse as well as in various astrophysical environments. In 
particular, HeH + was predicted to be observable in plan- 
etary and gaseous nebulae such as NGC 7027 0, 
However, despite searches for vibrational or rotational 
lines, no emission from HcH + has been detected from 
these objects so far, although it has been observed in 
laboratory plasmas for many years [H, Q . It was shown 
that the observation of the J — 1 — rotational line 
of HcH + is hindered by the near-coincidence of this line 
with a CH rotational line that has a greater intensity 
[Toj | . Other possibilities to detect HeH + include the po- 
tential presence of the ion in helium-rich white dwarfs 
[liT \v\ , metal-poor stars , or in supernovae [l4| , but 
the molecular ion has so far eluded observation. 

The formation of HeH + is mainly due to the radiative 
association between He and H + or between He + and H 
(processes {J) and (|5J)). While it has always been sup- 
posed that HeH + is formed in its ground X 1 S + state, 
one should also consider the possible role of the first 



metastable triplet state, a 3 £ + . This state, which corre- 
lates asymptotically to He + (ls) + H(ls), can indeed be 
populated and will not decay by collisions if the plasma 
density is low. As its radiative decay to the ground state 
is spin-forbidden, the state is metastable with a lifetime 
of 150 s [l5j. Moreover, in experimental studies on the 
dissociative recombination of HeH + [lj| [TtJ , it was sug- 
gested that the a 3 £ + state might be present in the ion 
beam and be responsible for a part of the cross section. 

In this work, we study the main mechanisms control- 
ling the abundance of HeH + in the a 3 S + state. The 
first reaction of interest is the photodissociation process, 
which destroys the molecular ion: 

HcH+(a 3 S+) +hu — ► (HcH+)* — ► He(lsnZ 3 L) + H+ 

— > He+(ls) + R(nl) 
(3) 

The photodissociation can occur following excitation to 
3 S + or 3 n electronic states and leads to atomic fragments 
He + H + or He + + H. We investigate this reaction with 
a time-dependent method based on the propagation of a 
wave packet on the excited electronic states coupled by 
non-adiabatic interactions that gives access to the con- 
tribution of each excited state to the cross section. 

Based on the results on the photodissociation process, 
we examine the radiative association of HeH + in the 
a S + state, which is one of the main mechanism of for- 
mation of this state. It occurs through a spontaneous 
radiative transition from the vibrational continuum of 
the excited b 3 S + state to a discrete vibrational level of 
the metastable a 3 S + state: 

He(ls2s 3 S) + H+ — ► HeH+(a 3 S+) + hv (4) 

The b 3 S + molecular state corresponds to collisions be- 
tween H + and He in its metastable ls2s 3 S state, which 
has a lifetime of about 8000 s [18( . Since metastable he- 
lium is easily produced in an ion source, the reaction (0| 
could lead to the presence of HeH + in the a 3 S + state 
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in the source. Since the b 3 £ + state has a lar ge p oten- 
tial well (0.73 eV), it could also be produced [19j . and 
its radiative decay would result in populating the a 3 E + 
state as well, a process that was investigated by Chibisov 
et al. [2(j. Metastable helium is also present in various 
astrophysical environments pH I22I ] , and in particular in 
planetary nebulae [23|, H3 , mostly due to the recombina- 
tion of He + with free electrons. The radiative association 
process ((4]) might therefore lead to the presence in space 
of HeH + in its triplet state. HcH + can also be formed in 
the a 3 £ + state by a transition from the continuum to a 
bound ro-vibrational state of the same potential energy 
curve. However, the rate coefficient for this process is 
extremely small [3]. 

We start in Sec. |n] by summarizing previous results 
of the molecular structure of HeH + as well as the theory 
of the photodissociation and radiative association pro- 
cesses. We present the photodissociation and radiative 
association cross sections and rate coefficients in Sec. IIII| 
and we discuss possible applications of the a 3 £ + state. 

II. THEORETICAL METHODS 
A. Molecular data 

We consider here all the 3 E + and 3 IT electronic states 
dissociating into atomic states with n = 1 — 3, where n is 
the largest principal quantum number of the atomic frag- 
ments. There is only one n = 1 3 E + state, while there are 
4 3 E+ and 2 3 IT n = 2 states, as well as 6 3 E+ and 4 3 n 
n = 3 states. This makes a total of 11 3 E + states and 6 
3 IT states. The adiabatic potential energy curves (PEC) 
for these triplet states are shown in Fig. [TJ The low- 
energy electronic spectrum of the molecular ion consists 
of states dissociating either into H + He + or into H + + 
He. For example, the lowest triplet state (a 3 S + ) disso- 
ciates into H(ls) + He + (ls), while the first excited state 
(b 3 S+) dissociates into H+ + He(ls2s 3 S). We adopted 
the molecular data presented in Ref. !25] . The potential 
energy curves, non-adiabatic couplings and dipole tran- 
sition moments were calculated at the complete active 
space self-consistent field (CASSCF) and configuration 
interaction (CI) levels. Since the excited electronic states 
undergo avoided crossings (see Fig. [TJ, the radial non- 
adiabatic couplings must be taken into account. Due to 
the large number of such couplings, we do not show them 
here and refer the reader to Refs. 0, [2(| where they are 
discussed. The matrix of the radial non-adiabatic cou- 
plings in the basis of the adiabatic electronic functions 
{Cm}, F mm > = (Cm\dn\C m i), can be used to build the 
diabatic representation [27|. The adiabatic-to-diabatic 
transformation matrix D(i?) is found by solving the dif- 
ferential matrix equation 

d R B(R) + ¥(R) ■ B(R) = (5) 

As in previous studies involving the excited states of 
HeH + |28l . |29| , we have only retained the couplings be- 
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FIG. 1. Adiabatic PEC of the n = 1-3 triplet states of HeH+ 
from Ref. [H]. Black, 3 E+ states; Red, 3 I1 states. The a 3 E+ 
and b 3 E + states are shown with their dissociation limits. 

tween adjacent states, F m . m ±i. Eq. ([5]) was solved nu- 
merically by imposing the initial condition ID(oo) = I. 
The diabatic potential energy curves are the diagonal el- 
ements of the matrix U d = D _1 • U • D, where U is the 
matrix of H el in the adiabatic representation, and the off- 
diagonal elements of U d are the diabatic couplings. On 
the other hand, we neglected the non-adiabatic rotational 
couplings, as their effect is expected to be negligible (30j . 



B. Photodissociation cross sections and rates 

The photodissociation cross section from a given rovi- 
brational state v, J is computed with a time-dependent 
wave packet method that has been described in detail in 
Refs. [28|, |29(, so we only briefly outline its main fea- 
tures. The first step is to calculate the rovibrational en- 
ergies E v j and wave functions ip v j(R) of the a 3 S + state, 
which is achieved using a B-spline method. The potential 
supports six vibrational states for J — 0, as was reported 
in Ref. [l5j]. The ground vibrational state is bound by 



3 



664.8 cm -1 while the v = 5 state is bound by less than 
2 cm -1 , and there is a total of 51 bound rovibrational 
states. 

The wave packet at t = 0, t — 0) is defined on 

each excited state / as the product of the wave function 
of the rovibrational state by the dipole matrix element 
fiif(R) = (Ci\d\Cf) between the initial state i and the final 
state /: 

& vJ {R,t = U) = l iif(R)il>vj{R) ■ (6) 

In this case, the initial state i is the a 3 S + state while 
the final state / can be any excited 3 £ + or 3 II state. 
The transition to £ states occurs for a parallel orien- 
tation of the field with respect to the molecular axis, 
while the transition to II states occurs for a perpendic- 
ular orientation. The wave packet subsequently evolves 
on the coupled potential energy curves, and the time- 
propagation is performed with the split-operator method 
in the diabatic representation [3l], [32[ . The total pho- 
todissociation cross section can then be computed from 
the Fourier transform of the autocorrelation function 

C{t) = E/<*L(^.°)I*L(^*)> P- The propagated 
wave packet therefore contains the information about the 
cross section for all energies. 

However, in this work we are not interested in the total 
cross section but rather in the partial photodissociation 
cross sections for each of the excited electronic states. 
We wish to understand which states form the dominant 
contribution to the cross section, and the partial cross 
sections are also necessary if one wants to compute the 
radiative association cross section (see below) . The com- 
putation of the partial cross sections is realized using 
a method based on the Fourier transform of the wave 
packet at a point Roq located in the asymptotic region 
[34| . The partial cross section for each final excited state 
/ is given by 

where /i is the reduced mass of the system, E p h = hv 

is the photon energy, kf — J 2^,{E v j + E p h — EL) is the 
wave number in the electronic state / with an asymptotic 
energy E[ s , and 

A* j = -j=J ^t j{Raojf)e ^ VJ+Epb )t dt (g) 

The calculations were performed using a spatial grid 
of 2 13 points from 0.5 to 100 au. The point at which the 
Fourier transform ([8} is evaluated was chosen as Roo — 75 
au. To avoid unphysical reflexions of the wave packet at 
the end of the grid, an absorbing potential starting at 
R = 80 au was introduced. The time step used in the 
propagation was 1 au, and tests with shorter steps were 
performed to assess the convergence of the cross sections. 
For the highest vibrational level, v = 5, we extended the 
grid to 150 au and took R^ = 100 au. This is necessary 



as the wave function for this weakly bound level extends 
to large internuclear distances. 

The rate coefficient for the photodissociation following 
the absorption of a photon emitted by a central star of 
radius i?* with blackbody temperature at a radial 
distance R is given by Q 

fc ( T *) = ^ / e E pl f k p B T> n ^pO dEph ' 

(9) 

We have chosen the value (R+/R) 2 = 10 13 , as suggested 
in Ref. . 



C. Radiative association 

The cross section for the radiative association process 
can be obtained from the photodissociation cross section 
using the detailed balance principle. The radiative asso- 
ciation cross section for the formation of the a 3 E + state 
into a level v, J alon g th e potential of any excited state 
can be calculated as [35j 

<j(Ek)^-^r^j(E pil ) . (10) 

In this equation, &yj(E p h) is the partial cross section 
corresponding to the photodissociation of the v, J level 
of the a 3 S + state via an excited state as a function of 
the photon energy, while cr^j{Ek) is the cross section for 
radiative association into the v, J level as a function of 
the relative kinetic energy of the incident particles. It is 
related to the photon energy by 

E k = E ph - AE , (11) 

where AE is the difference between the dissociation en- 
ergy of the electronic state along which the collision takes 
place and the energy of the rovibrational level of the 
a 3 E + state in which the association is realized. As 
the association can occur into any bound rovibrational 
level, the cross section for radiative association must be 
summed over all possible values of v and J. It should 
be noted that, by construction, the radiative association 
cross section from a particular electronic state calculated 
using Eq. (fTTU) includes non-adiabatic effects. 

The radiative association rate constant is calculated 
assuming a Maxwell-Boltzmann distribution of incident 
particles, 

k(T) = (r^) 3/2 477 r E k e- E ^ T a{E k ) dE k . 
\k B TJ ^AmJo 

(12) 
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III. RESULTS AND DISCUSSION 
A. Photodissociation 

The photodissociation cross section for the initial rovi- 
brational level v — 0, J = is presented in Fig. O The 
total cross section, as well as the contributions of the 
3 £ + and 3 II states, are shown. In the case of dissocia- 
tion through 3 S + states, the cross section is composed 
of two peaks centered around 9.5 eV and 11.6 eV. The 
first peak is due to the b state and the second to the 
three remaining n = 2 states. The small structure in the 
range 15 — 17 eV is due to the n = 3 states. We note 
that the cross section starts abruptly at 8.92 eV, which 
corresponds to the energy difference between the initial 
state (i.e., the v = level of the a 3 S + state) and the 
dissociation energy of the b 3 S+ state. The b 3 S+ state 
has a deep potential well of about 0.73 eV, as can be seen 
from Fig. [T] with an equilibrium position R e = 7.75ao 
125] • A component of the initial wave packet © will have 
an energy below the dissociation energy of the b 3 E + 
state, and this component will oscillate in the potential 
well. Therefore, it never reaches the asymptotic region 
and makes no contribution to the cross section of the 
b 3 S + state. On the other hand, the component of the 
wave packet with sufficient energy will travel toward the 
asymptotic region, which explains the threshold behavior 
of the cross section. As in the case of the photodissocia- 
tion from the ground X 1 T, + state of HeH + (28|, we ob- 
served a strong influence of the non-adiabatic couplings 
on the photodissociation cross sections, due to the large 
number of avoided crossings affecting the excited states 
PEC (see Fig. [TJ. 

The cross section for dissociation via 3 II states is dom- 
inated by a narrow peak centered at 10.3 eV and a wider 
peak around 13 eV, both due to the two n — 2 states. 
The shape of the cross section can be explained using 
the reflection principle, which states that the cross sec- 
tion reflects the spatial distribution of the wave function 
of the initial state, and that its width is proportional to 
the steepness of the potential energy curves in the ex- 
citation region 36]. As can be seen in Fig. [TJ in the 
excitation region the PEC of the lowest 3 n state is al- 
most flat, leading to a very narrow cross section, while 
that of the second 3 n state is steeper, resulting in a wider 
cross section. The diffuse structure in the cross section 
between 14 and 17 eV is due to the n = 3 states. As these 
states interact strongly through non-adiabatic couplings, 
the cross section is much less structured than for the two 
n = 2 states. 

The photodissociation cross section is strongly depen- 
dent on the initial vibrational level. This is illustrated in 
Fig. [3J which shows the cross section for the initial state 
v = 1, J = 0. It has a more complicated structure than 
for the v — case, reflecting the oscillatory behavior of 
the initial wave function, but we still observe the domi- 
nance of the 3 n states as well as a threshold behavior at 
an energy corresponding to the dissociation energy of the 
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FIG. 2. Photodissociation cross section for the initial state 
v — 0, J = 0. Blue dashed lines: contribution of the 3 E + 
states. Red dotted lines: contribution of the 3 II states. Black 
full line: total cross section. 
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FIG. 3. Photodissociation cross section for the initial state 
v — 1, J = 0. Blue dashed lines: contribution of the 3 E + 
states. Red dotted lines: contribution of the 3 II states. Black 
full line: total cross section. 



b T, + state. We reach similar conclusions for the higher 
vibrational levels v > 2. On the other hand, the cross 
section is almost independent of J and is insensitive to 
rotational excitation as the initial wave function depends 
only weakly on J. 

The photodissociation rate coefficient, given by Eq. 
©, is presented in Tab. [J as a function of the star tem- 
perature T* and for three different value of the matter 
temperature T m . The matter temperature T m governs 
the rovibrational distribution of the initial state of HcH + . 
Assuming equilibrium, the population of each state is 
given by 

PvJ {T) = (2 J + 1) exp{-E vJ /k B T m )/Z(T m ) (13) 
where ks is the the Boltzmann constant and Z(T m ) is 
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the partition function. The photodissociation cross sec- 
tions and rate coefficients must take this dependence into 
account. We give in Tab. Q]the contribution of the b 3 E + 
state to the rate coefficient (denoted by kb), as well as 
the total contribution of the 3 £ + and 3 II states (ks and 
kn, respectively) and the total rate (fctot)- We observe 
that even at low temperature, it is not sufficient to con- 
sider only the b 3 E + state to obtain an accurate value 
of the photodissociation rate coefficient. Other excited 
electronic states, and the 3 I1 states in particular, must 
also be taken into account. For example, at T+ = 10 4 K 
the contribution of the b 3 S + state to the total rate co- 
efficient is less than 50% while at T+ = 5 x 10 4 K it does 
not exceed 20%. This is in stark contrast to the pho- 
todissociation from the ground state of HeH + , in which 
case the rate is dominated by dissociation into the first 
excited state (37j | . 



B. Radiative association 

The radiative association cross section along the b 3 S + 
state, calculated using Eq. (fTU)) . is presented in Fig. @] 
as a function of the collision energy. The cross section is 
dominated by transitions to the levels with v = — 2, and 
it is several orders of magnitude larger than the cross sec- 
tion for radiative association along the a 3 S + state that 
was studied Kraemer et al. [4J. In that case, the as- 
sociation takes place via a radiative transition from the 
continuum to a bound level of the same potential, lead- 
ing to a small cross section. The cross section shown 
in Fig. [4] decreases sharply for energies over 2.77 eV, 
which corresponds to the threshold energy above which 
the transition probability from the b 3 S + state to the 
a 3 S + state becomes strongly suppressed. It should be 
noted that the cross section presented in Fig. |4]is only the 
non-resonant part of the cross section and that the contri- 
bution from shape resonances due to the centrifugal po- 
tential, as well as resonances caused by the non-adiabatic 
interactions, should also be included. Unfortunately, the 
accurate determination of these resonances requires im- 
practically long propagation times [19j . Taking shape 
resonances into account results in an enhancement of the 
radiative association rate coefficient and can be evaluated 
in a separate calculation using the Breit-Wigner theory 
[HI, Hll . The contribution of resonances on the radiative 
association rate can be quite large in cold environments, 
but it decreases quickly with increasing temperature [4(| ■ 
In addition, since the a 3 S + state has a small binding en- 
ergy and supports only 51 bound rovibrational states, we 
can reasonably expect the number of shape resonances to 
be small and their contribution to the rate coefficient to 
be negligible except at very low temperature, similarly 
to what was recently shown for the radiative association 
of LiHe+ gjj. 

The rate coefficient, calculated using expression (IT21 . 
is shown in Fig. [SJ It presents a maximum at 10 K 
and decreases at higher temperature. It is smaller than 
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FIG. 4. Radiative association cross section for the radiative 
association along the b 3 E + state, process (TjTJ). 
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FIG. 5. Rate constant for the radiative association process 
b 3 E+ -> a 3 E+ + hu. 



the rate coefficient for association in the ground state 
However, the value of the rate is probably slightly 
larger than the one presented in this work due to the 
contribution of shape resonances that we neglected. 

As we have calculated the partial photodissociation 
cross sections into all the n = 1 — 3 excited states, we 
can easily evaluate the cross section for radiative asso- 
ciation in a discrete level of the a 3 S + state following 
a spontaneous transition from the continuum of any of 
the excited electronic states. In Fig. [51 we show the 
rate coefficients for the dominant channels of formation 
of HeH + in the a 3 E + state. It should be noted that 
it is necessary to include the effect of the non-adiabatic 
couplings to obtain accurate rate coefficients, due to the 
strong interaction between the excited electronic states 
of HeH + . In the 3 E + symmetry, the radiative associa- 
tion process is dominated by formation along the b 3 S + 
state. The electronic channels that give the next largest 
rate coefficients for the radiative association of the a 3 S + 
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TABLE I. Photodissociation rate coefficients (s -1 ). T m (K) is the matter temperature and 7*(K) is the temperature of the star. 
kb denotes the contribution of the b 3 E + state to the rate coefficient, fcs is the contribution of all the 3 E + states, kn is the 
contribution of the 3 fl states, and k to t — ks + kji is the total rate. 
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state are the second and fourth n = 2 states (see Fig. 
[T]), which correlate asymptotically to He(ls2p 3 P°) + H + 
and Hc + (ls) + H(2s), respectively. However, the rates 
are much smaller for these states. On the other hand, the 
approach along the two lowest 3 LT states (dissociating into 
He(ls2p 3 P°) + H+ and He+(ls) + H(2p), respectively) 
yields rate coefficients that are much larger than the rate 
corresponding to the transition b 3 S + —> a 3 £ + . This 
shows that radiative association along the b 3 £ + is not 
the most efficient way of producing HcH + in the a 3 E + 
state. 



C. Applications 

Metastable helium can be present in significant con- 
centration in astrophysical environments such as plan- 
etary nebulae due to recombination of He + with elec- 
trons (42j, and it is conceivable that it might influence 
the abundances of various atomic and molecular species. 
Roberge and Dalgarno [2| already considered the role of 
metastable helium on the production of HeH + in the 
ground state following collisions between He(ls2s 3 S) 
and atomic or molecular hydrogen and showed that under 
some circumstances, these reactions can be a substantial 
source of HcH + . 

Based on the rates presented in Sec. IIII1 we can pro- 
vide a rough estimate of the abundance of HeH + (a 3 E + ) 
in planetary nebulae. If we assume equilibrium be- 
tween photodissociation and radiative association along 
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FIG. 6. Rate constants for radiative association of HeH in 
the a 3 E + state along various electronic states. Black full 
line: association along the b 3 E + state; dotted blue line: as- 
sociation along other 3 E + state; dashed red line: association 
along 3 fl states. 

the b 3 S + state, we have 

n(HeH+)fc d = n(2 3 S>(H + )fc a (14) 

where n(HeH + ), n(2 3 S), and n(H + ) are the density of 
HeH + , He(ls2s 3 S') and H + , respectively, while fc d and 
k a are the rates for the photodissociation and radiative 
association processes. We consider here a typical neb- 
ula with parameters n(H + ) = 10 4 cm" 3 , n(He + ) = 10 3 
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cm" 3 , T* = 5 x 10 4 K and T e = 10 4 K [23[. A formula 
giving the abundance of metastable helium in nebulae as 
a function of the He + density and electron density and 
temperature was derived by Clegg (24{. Using the pa- 
rameters above, we get a density of metastable helium 
n(2 3 S) ~ 4 x 10~ 3 cm -3 . Combining with the results 
for the rates fc d and fc a given respectively in Tab. U and 
Fig. G3 we get n(HeH+) = 3 x 10~ 10 cm" 3 . The den- 
sity of HeH + in its triplet state is much smaller than 
in the ground state but it can nonetheless influence 
the abundance of other species through reactions ([3]) and 
(U]). The radiative association process ((H) could also be 
of importance in the chemistry of the early universe [![. 
Neutral helium was first produced following recombina- 
tion of He + , therefore populating the 2 3 S state. While 
metastable helium was probably re-ionized due to its low 
binding energy, it could still form HeH + by radiative as- 
sociation with H + . 



IV. CONCLUSIONS 
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its metastable triplet state by means of time-dependent 
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